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We propose the development of an instrument by the Martin & Puplett-type 
Fourier Transform Spectrometer to applying the aperture synthesis technique in 
millimeter and submillimeter waves. We call this equipment the Multi-Fourier 
Transform interferometer(MuFT). MuFT performs a wide band imaging, spec¬ 
troscopy and polarimetry in millimeter and submillimeter wavelengths. We describe 
the fundamentals of MuFT, and give an example of one potential implementation. 
Full description of the observables by MuFT are provided. A physical explanation 
of the observability of the complex visibility by MuFT is given. Fundamental re¬ 
strictions on observations with MuFT, eg. limits on spectral and spatial resolutions 
and field-of-view, are discussed. The advantages of MuFT are also summarized. © 
2008 Optical Society of America 
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1. Introduction 

Astronomical observation of submillimeter waves represents one of the last unexplored regions of 
the electromagnetic spectrum. One of the difficulties with performing observations in this band is 
owing to a large atmospheric absorption. Since submillimeter waves are a boundary of the radio and 
infrared, there remain a lot of frontiers in the development of fundamental observational technology. 

Radio interferometry and synthesis arrays, which are basically ensembles of two elements inter¬ 
ferometers, are used to make measurements of fine angular resolution. 1 Interferometers have been 

1 When this research was performed, Izumi S. Ohta was with Tohoku University Aoba Aramaki Aoba-ku Miyagi 
980-8488,Japan. Izumi S. Ohta is now with National Astronomical Observatory 2-21-1 Osawa Mitaka Tokyo 181-8588, 
Japan 
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used at millimeter waves 2 and have been considered one of the best instruments at submillimeter 
wave. 3 By setting two telescopes with separation larger than the maximum construct-able diameter 
of the single mirror, higher angular resolution than single dish system is easily achieved with aper¬ 
ture synthesis type interferometers. By combining many telescopes, a large effective area can be 
attained to increase the sensitivity of the system. Interferometers are robust against atmospheric 
fluctuations since large scale atmospheric fluctuations do not cause interference signals and are 
automatically omitted by the interferometers. 

However, currently operating millimeter and submillimeter interferometers have the following two 
fundamental problems because they use heterodyne receivers as focal plane detectors: The receiving 
systems perform modulation and amplification of the source signal by mixing with a local oscillator 
signal. They can measure the phase of the incident electromagnetic wave. In the following discussion, 
the space-borne observation assumes as for an ideal case. The noise temperature of heterodyne 
receivers at frequency v is written as Tsys-het ~ 2(hu/k B ) \(e 2 n 2 (v,T 0 ) + en(V, T 0 )) 2 + ll, where n 
is the thermal Planck function with To = 2.73 K due to the cosmic microwavebackground (CMB), 
h and k B are Planck constant and Boltzmann constant, respectively, and e represents the detection 
efficiencies. The first term represents the photon noise due to statistical fluctuation of the CMB 
intensity. The second term represents the quantum limit of the heterodyne receivers which comes 
from the uncertainty principle between the phase and the number of photons. 4 This equation 
shows that the noise temperature is limited by the quantum limit and is increased linearly with the 
frequency in the frequency range of v > 100 GHz. Therefore, the sensitivity of heterodyne-based 
interferometers is severely limited by the quantum limit of the receivers. Another problem is the 
severely limited field-of-view (FOV) of the interferometers, due to the difficulty of constructing a 
large format detector array of heterodyne receivers. 4 In addition, the possible bandwidth of the 
interferometer is also limited by a bandwidth of Intermediate Frequency(IF) amplifier, since the 
phase change of the modulated signal after mixing with a local oscillator signal must be followed. 

On the other hand, direct detectors such as bolometer are also used in millimeter and submil¬ 
limeter systems except for interferometers. Direct detectors measure the intensity of the source 
signal. They cannot measure the phase of the incident electromagnetic wave. Therefore, they have 
not been used as focal plane detectors for interferometers. The quantum limit of a direct detector 
is the minimum detectable power when any internal noise can be neglected, the quantum limit of 
a direct detector is expressed in noise temperature as T qu _di r ~ hv j(k B AtAu), where At and Az/ 
are the response time of the detector, which is about 1 — 100 msec in the case of bolometer, and 
the bandwidth of the detector, which is about 10 — 100 GHz. 4 The noise temperature of direct 
detectors at frequency v is written as Tgys-dir ~ 2 {hv/k B ) |(e 2 n 2 (z/,T 0 ) + en(v,T 0 )) 2 + 1/AtAz/] 
The scecond term is negligible compared with the CMB photon noise at millimeter and submi 
limeter waves. Therefore, noise temperature of direct detectors are much lower than heterodyne 
receivers in the frequency range of u > 100 GHz. Actually, the noise power of the most sensitive 
currently operating bolometers is already lower than the quantum limit of the heterodyne in the 
frequency range of v > 100 GHz. In Contrast to the heterodyne, a large format millimeter and 
submillimeter direct detector array is possible to construct. 5,6 There is essentially no limit on the 
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bandwidth for direct detectors. 7 

The above discussion indicates that the bolometric interferometer, which is able to use a direct 
detector as a focal plane detector for the interferometer, could be one of the best ideal instruments 
for millimeter and submillimeter astronomy since it shares the advantages of both direct detectors 
and interferometer. Since direct detectors cannot obtain the information of the phase of electromag¬ 
netic waves, two beams captured by two apertures must be mixed before guided into the detectors. 
There are two ways to mix beams. One is called Fizeau type beam mixing. The simplest example 
of this method is that two beams obtained by masking a single telescope are mixed at the focus of 
the telescope. Images of different portions of the sky are focused on different positions of the focal 
plane. Therefore, the extension of the FOV by using the focal plane detector array is straightforward 
in this case. However, there is a fundamental problem in applying this method to millimeter and 
submillimeter observations. We have to use many detector pixels within a single FOV to resolve the 
fringe pattern inside of the FOV. Since the direct detector arrays in millimeter and submillimeter 
wave bands are still expensive, the cost performance of applying this method to this waveband is 
poor. The other method is called the Michelson-type beam mixing. An application of the Fourier 
Transform Spectrometer (FTS) to aperture synthesis interferometry has been studied as a possible 
solution of the Michelson-type beam mixing bolometric interferometer. The application was first 
independently proposed by Itoh & Ohtsuka(1986) 8 and Mariotti & Ridway (1988) 9 in near infrared 
(NIR). The system was referred by these authors as ’’double Fourier interferometry”, since the 
signal obtained by the system is a Fourier transformation of both spectra and intensity distribu¬ 
tion on the sky. Itoh & Ohtsuka studied single pupil interferometry. Mariotti & Ridway studied 
multi-pupil interferometry for high spatial resolution. As explained in this paper, the extension of 
the FOV by putting the focal plane detector array is not straightforward, but it is possible with 
some restrictions. Rinehart et al. (2004) 10 have been studying the extension of the FOV by putting 
a detector array on the focal plane of the double Fourier in laboratory experiments in an optical 
wave bands by using optical CCD camera. They refer to this experiment as the Wide-held Imaging 
Interferometry Testbed. 

In this paper, we propose the application of a Martin & Puplett-type Fourier Transform spectrom¬ 
eter (Martin & Puplett(1969), 11 hereafter MP-FT) to the aperture synthesis system in millimeter 
and submillimeter waves(0hta(2004) 12 ). A wire grid polarizer (WG) is used as a beam splitter in 
this system. So, the wavelength dependence of reflectivity of WG is small and is suitable for a wide 
band measurement system. By setting the two input WGs appropriately, 2D intensity distribution 
of four Stokes parameters can be measured. The signal obtained by our system is a Fourier trans¬ 
formation of spectra and intensity distributions of four Stokes parameters (multiple components) 
on the sky. Therefore, we refer to this system as a Multi-Fourier Transform interferometer, abbrevi¬ 
ated to MuFT. The abbreviation also contains the meaning that this instrument measures Mutual 
correlation of the source signal instead of auto-correlation as in the case of usual FTS. 

The plan of this paper is as follows: In Section 2, fundamentals of imaging and spectroscopy by 
MuFT are explained using scalar waves for simplicity. How complex visibility is measured by MuFT 
is also explained in Sec. 2. In Sec. 3, details of the components of MuFT are introduced. A full 
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description of the observables by the MuFT and restrictions intrinsic to MuFT are also provided 
in Sec. 3. In Sec. 4, the advantages of the MuFT are summarized. A summary of the paper is given 
in Sec. 5. This paper presents the fundamental theory of MuFT. Details of broadband imaging and 
spectroscopy experiments using MuFT are presented in a forthcoming paper. 


2. Fundamentals of Multi-Fourier Transform interferometer 


2. A. Fundamentals of imaging by MuFT 

Coordinate systems on the source plane and the observer plane are defined in Fig. 1. The notation 
of this coordinate system follows Born & Wolf . 13 The coordinate system on the source plane (£, rf) is 
fixed on the sky. The origin of the coordinate system on the source plane is defined to the center of 
the FOV. The coordinate system on the observer plane is defined as X and Y axes that are parallel 
to £ and 77 axes, respectively. The origin of the observer coordinate is defined as the mid-point of 
two apertures mounted at Pi and P 2 . The position vector connecting to the center of the aperture 
mounted at P 2 with the center of the aperture mounted at Pi is defined as b = (b x ,b y ), which 
is referred to as the baseline vector. No absorption, no emission and no scattering due to matter 
distributed between the source and observer are assumed. 

We consider the measurement of the mutual coherence function of the waves obtained by two 
apertures Pi and P 2 . The source electric fields arrived at Pi and P 2 are defined as Ep\(t) and 
Ep 2 (t), respectively. They are recombined at Q. When two beams are recombined at Q, there is 
a time lag r which is caused by the internal light path length difference from each aperture to Q, 
that is r = (Jt\ — £ 2 )/c = 2x/c. The mutual coherence function with the time lag r is defined as 

Fi 2 (t) = {Epi(t)Ep 2 (t + t)) , (1) 


where (} denotes time average. The intensity of the recombined signal is written as 


Ii2(Q,r) = /i+/ 2 + 2rS r 2 ) (r), 


( 2 ) 


where I\ and I 2 are the intensities of the signals obtained by each aperture and F^ (t) is a real part 
of ri 2 (r). The unwanted DC components could be subtracted, and are neglected in the following 
discussion. Since the sizes of the astronomical objects are much larger than the wavelengths of the 
electromagnetic waves which we are interested in, any astronomical extended source can be treated 
as a collection of incoherent point sources except for a few exceptional cases. A viewer angle of 
the source from the observer plane is introduced as 0 = (0 x ,9 y ) = (£/R,rj/R). Then, the mutual 
coherence function and its real part can be expressed by the specific source intensity, 1(0,u), as 


Fi 2 (b,r) = ri 2 (r) = 1(0,v) exp 


— 2irj — (b • 0 ) + 27 xjvr 


rSJtb ,r) = r 


= — 


12 


T = 


1(0 , v) cos 


— 27T — (b • 0 ) + 27 TUT 


dnjd z 6, 
diA d 2 e, 


( 3 ) 

( 4 ) 


where D is the solid angle expanded by the source, j is an imaginary unit, and T^b, r) and 
rg(b, r) explicitly express the dependence of the mutual coherence function on the baseline vector. 
The intensity expressed in Eq. (4) is the observable by MuFT. 
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(r) 

The spectroscopy-resolved source image is extracted from the observed r^ 2 in the following way: 

(r) 

For each fixed combination of b, T^ 2 is measured for the time lag from —To to to by changing 1 2 

(r) 

continuously. By performing a Fourier integration of the obtained T^ 2 in r, the spectrally resolved 
mutual coherence function, is), is obtained as 


12 


(h,is) = J {/ 1(6, is')2tqS\uc \2ir(is' — is)tq\ exp 


—27rj — (b • 6 ) dis' 


d 2 e, 


( 5 ) 


where 1(6,—is) = 1(6, is) is used and sine 2 ; = (sinx)/x. This is a kind of Winner-Khenchine 
Formula. 13 The integration by is 1 reduces to is' = is with spectral resolution of Ais ~ 1/ro- Since 
is 1/ro in an usual case, Ais <c is. Now the customary notations in the interferometer can be 
introduced, that is 


b x v 


VE" u X bylS Vy 

u = -= - r ,u= —= -f. (6 

c A c A 

The (u,v) is called spatial frequency. When |b • 6\Ais/c -C 1, the Fourier transformation of the 
mutual coherence function can be written as 


fu(u,v,v) = J^i(6,v)exp[-2irj(u9 x + v6 y )] d 2 0, (7) 

where I is a spectrally convolved source intensity distribution defined as 

1(6,is) = /l( 0 ,, / )2rosinc[27r(,'-,)To]^. (8) 

The functional form of the convolution kernel depends on the apodization function adopted when 
the Fourier transformation by r is performed. In the above example, the top hat function is adopted 
as the apodization function. Eq. (7) is equivalent to the equation which expresses the van Cittert- 
Zernike Formula in Fraunhofer limit. 13 It shows that the Fourier transformation of the mutual 
coherence function, T\2(u,v,is), is the complex visibility function of the source except for the 
proportionality constant. The main reason why the MuFT is capable of measuring the source 
complex visibility function, is to take data about the time lag r from —r D to r D symmetrically. 
By Fourier transforming them in u and v, source images are obtained for various frequencies. The 
MuFT is, in principle, is capable of performing simultaneous measurements of source image and 
spectrum. The MuFT corresponds to the XF-type radio interferometer. 1 Physical explanation of 
the principles of imaging by MuFT is given in Sec.2B. 

The observability of the complex visibility function by MuFT makes observations with MuFT 
easier. From Eq. (7) T\ 2 (u, v ,v) = Ti 2 (—u,—v,is). Here we have used the fact that the source 
intensity distributions, I (6, is), are a real variable called the reality condition. Therefore, sampling 
in only half of the uv plane is enough to perform imaging observations by MuFT. 


Large dynamic range advantage When the source intensity distribution does not depend 
on frequency, all the source images obtained for various frequencies can be summed up into a 
single image in the following way: In this case, the frequency dependence i(is) and the spatial 
coordinate dependence of the source intensity B(6) are decomposable, such as 1(6, is) = B(6)i(is). 
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The spectrally convolved source intensity distribution can also be written as I(6,v) = B(9)i{v) 
where i(y) is the convolution of i(y) with the spectral convolution kernel. Eq. (7) can be rewritten 
as 


ti 2 (u,v,v) = / B(6) i(i/)exp 


-2nj-(b ■ 0) 


d 2 o. 


( 9 ) 


Divided by z(z/), the frequency dependence of the source intensity can be removed. 

Suppose that the observations for the baseline vector combinations on a circle with a radius of 
|b| = b on the baseline vector plane were performed with continuous sampling. Observed frequency 
bandwidth is supposed to be very wide. Adopt a polar coordinate on the baseline vector plane as 
b = (6, ip). Integrate T^b, v) exp[27rjV(b-0)/c]/i(V) by vdv and ( b 2 /c 2 )dp> for a fixed b. Transform 
variables from v and ip to u = bvcostp/c, v = busiiup/c using d(u,v)/d(v,<p) = vb 2 /<?. This 
procedure leads to the source intensity distribution as 


[ 2n \ r ri 2 (b,i/) 

/ 1 / v -■( \ exp 

/ o Jo l(v) 


2 vrj-(b • 6) 


du \ ^-dp 


r c roc 

I o [Jo 


V f d 2 9'B(0 / ) exp 

Jn 


exp 


2 vrj-(b • G) 


—2irj— (b • O' 

J (/ B(O') exp ^2irj u(Q' x — 9 X ) + v{0' y — 9 y ) | d 2 ^ dudv, 

[ B(0')5 2 {6'-6) d 2 9', 

J n 

B(G). 




( 10 ) 


In essence, the procedure can be thought of as follows: Since the spatial distribution is independent 
of frequency, the choice of frequency to measure the source intensity distribution is arbitrary. For 
a fixed baseline interval, observations with short wavelengths measure precise structures with high 
resolution and observations with long wavelengths measure overall structures with low resolution. 
By gathering all this information, precise source intensity distribution is obtained. 

Several modifications are required to apply the above procedure to real observations. First of all, 
sampling on the baseline vector plane can neither be perfect nor continuous. The sampled points 
on this plane will be discrete and sparse. Instead of the Fourier integration described in Eq. (10), 
discrete Fourier transformation must be performed. To express the sparseness of the sampled points, 
convolution with a window function which is set at zero at the baseline vector points where the 
observation was not performed, must be taken into account in this calculation. Secondly, the source 
spectrum may not be finite for all frequencies. The source intensity may be zero in some frequency 
ranges due to absorption by the intervening medium or no emission from the source. Division by 
i(v) described in Eq. (10) will diverge in these frequency ranges. The data for these frequency ranges 
should be removed in the procedure. Otherwise, they may introduce a large amount of noise into the 
obtained map. These modifications result in an imperfect map. However, it does not mean that the 
intrinsic power of the large dynamic range advantage is lost. The imperfection of the synthesized 
map due to the imperfect uv coverage is one of the most common problems in interferometer. The 
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large dynamic range advantage of MuFT can be explained as a large uv coverage being obtained 
from a small number of baseline vector combinations. This technique is essentially the same as the 
Multi-Frequency Synthesis (MFS) used in radio interferometers (Conway et al.(1990) 14 ). 

The situation to which this procedure is applicable is a very special case. However, there are plenty 
of astronomically important objects whose spatial distribution is independent on frequency, such as 
millimeter and submillimeter wave intensity distribution of clusters of galaxies known as Sunyaev- 
Zel’dovich effect (hereafter SZ effect, Sunyaev &; Zel’dovich (1972), 15 Hattori & Okabe(2005) 16 ). 

2.B. The physical explanation of imaging by MuFT 

As shown in the previous section, the essential reason why MuFT is capable of performing imaging 
observations in broadband is that the complex visibility function of the source can be extracted 
from the observables with MuFT. This is accomplished by changing the time lag from — To to to 
symmetrically. The reason why the complex visibility function is measurable by this process, is 
physically explained as follows. 

Consider a broadband point source observation by MuFT. Suppose that there is a geometrical 
time lag, r g , caused by the light path length difference from the source to each aperture Pi and P 2 , 
that is T g = (R\ — R 2 )/c. A burst-like interference signal is observed as scanning in r by changing 
£2 continuously at r = — r g , where the geometrical time lag is compensated by the internal time 
lag. Since this burst-like signal always appears at r = 0 in the FTS, the signal appearing at r = 0 
is called a central burst of the interferogram. We also refer to this burst-like signal observed for 
a point source in MuFT to a central burst as following the FTS. In the case of MuFT, the burst 
center appears at the position where the total time lag of the light led by each aperture is zero. 
The position of the central burst observed by MuFT contains information about the position of the 
point source in the sky. If the source position is slightly shifted from the origin of the source plane 
toward the same side of the aperture Pi and r g is negative, then r must be positive and vice versa 
to observe the central burst signature. Suppose that r is only scanned in the positive range. The 
central burst of the second source is not observable by this scan and we cannot know the position 
of the second source. By rotating the baseline vector by 180 degrees if possible, the geometrical 
time lag for the second source becomes negative and the central burst caused by the second source 
becomes observable by the positive r scan. On the other hand, the central burst of the first source 
goes out of the positive r region. This rotation corresponds to take a —b configuration of two 
apertures. Therefore, when r is scanned only in the positive range, to obtain a full source image, 
measurements with b and — b configurations are required. However, if we allow for scanning of the 
negative r region, the central burst of the second source can be detected in the negative r region 
without rotating the baseline vector. Therefore, by performing a scan in r from the negative to the 
positive range symmetrically, the full source image is measurable without a performing 180-degree 
rotation of the baseline vector. This is why the symmetrical scan in t leads us to obtain the complex 
visibility function. 

In contrast to the MuFT, the complex visibility function is obtained in the following ways in 
the heterodyne system: The real part of the complex visibility function is obtained by directly 
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correlating the signals obtained by each aperture. The imaginary part is obtained by correlating 
the signals after applying the 7r/2 phase shifter for one of the signals. 

3. The application of the Martin Puplett-type FTS to aperture synthesis 

3.A. Components of MuFT 

Fig.2 is a schematic diagram of the MuFT system. From the top of the page, the parts depicted 
are the Light Collecting unit(LiC), the Fourier Interference unit(FI) and the Detection & Sampling 
unit(DeS). The detected signal is transfered to a computer after it has been through an electrical 
amplifier, a noise filter and an AD converter. The sampled data are analyzed by a data analysis 
system. Since the data is transformed into Polarimetry, Spectroscopy and Imaging information in 
the data analysis system, the data analysis system is referred to as PSI. 

The LiC first part performs a division of the wave front of the incident source signal, and guides 
the signal to the Fourier Interference unit. Two different types are considered. One type has multi¬ 
apertures installed on separated multi-platforms. The antenna in multi-platforms type radio inter¬ 
ferometer and siderostat in multi-platforms type optical-IR interferometer correspond to this type. 
In this type, there are large geometrical time lags except for the moment when the source is at 
southing. To compensate for the geometrical time lags which are also time-variable due to earth ro¬ 
tation, a large and time-variable offset in the scanning center of r must be introduced in MuFT. To 
do this, a fringe tracking delay-line technique can be applied. 1 ' When the MuFT system is applied 
to multi-platforms type LiC, performing a symmetrical scan in r from the negative to the positive 
region is essential to obtain the complex visibility function. This is because the possible rotation 
angle of the baseline vector relative to the astronomical sources is limited, since the sources set and 
are not observable for nearly half a day except for sources around a celestial pole. The other type 
of LiC has that multi-apertures are installed on a single platform. For simplicity, we consider the 
two-aperture case. The positions of the two telescopes are controlled by controlling the azimuth and 
elevation angles of the platform so that a fixed reference point on the targeted source is maintained 
by a baseline vector bisector line, perpendicular to the baseline vector. The same reference point 
must be kept in the center of the FOV of each telescope. Examples of this type of LiC are CBI 18 
and AMiBA 19 which have multi apertures installed on a single platform. Two mirrors connected 
by one rail mounted on a turntable, such as SPIRIT (Mather & Leisawitz(2000) 2 °), could be a 
solution. In the space, using multi-satellites could be an attractive solution. 

The Fourier Interference unit is an instrument in which divided light beams are recombined after 
modulating one of the light path lengths to make an internal time lag and is the heart of the MuFT 
system. The method of beam mixing is the Michelson type, as mentioned in the introduction. This 
is designed by applying the MP-FT to aperture synthesis. It has two entrance windows to receive 
two light beams guided from the light collecting unit. Wire grid polarizers are used as beam splitters 
and beam combiners. One example of the FI unit is shown in Fig.3. 

The direction of the WGs 1 and 2 wires must be set in the following ways. Consider light rays 
which enter the WGs 1 and 2 along the vertical axis. These rays are reflected by the WGs and travel 


toward the Roof Top Mirrors (RTM) along a horizontal axis, and are then reflected by RTMs and 
travel back toward the WGs 1 and 2. Consider the projection of the WGs 1 and 2 on the planes 
normal to these light rays. The wires of these WGs must be tilted 45 degrees from the vertical axis 
on these planes. There are two possibilities for the wires being tilted clockwise (CW) or counter 
clockwise (CCW) from the vertical axis. There are four combinations for the alignment of the WGs 
1 and 2. Two basic configurations of the WGs are illustrated in Fig.4. When the wires of both 
WGs are set in a CW direction, Stokes parameters of X and Q are observable. This is referred to 
as option 1. By setting both WGs to CCW, complementary data for the first case is obtained. One 
of the configurations gives X + Q and the other gives X — Q. When the wires of WG 1 are set CW 
and those of WG 2 are set CCW direction, Stokes parameters of U and V are observable. This is 
referred to as option 2. By changing the setting of WGs 1 and 2 vice versa, complementary data for 
the original combination are obtained. One of the configurations gives U + V and the other gives 
U — V. Therefore, by performing observations with all four combinations, spectrally resolved source 
intensity distributions of 4 Stokes parameters are obtained. Since only three of 4 Stokes parameters 
are independent for completely polarized waves, systematic errors could be self calibrated in MuFT 
by measuring all 4 Stokes parameters for the completely polarized light beams. For the intensity 
measurement, observations with option 1 is enough. Details are explained in the next subsection. 

The other WGs installed in the FI unit for each option must be aligned in the following ways: 
In option 1, the wires of BS WG must be directed vertically or horizontally as illustrated in Fig.3 
to combine the two beams taken by WGs 1 and 2 at BS. The combined beam, either reflected or 
transmitted by BS, is led to WGs A or B. To get interference terms of the two beams, the wires 
of output WGs A and B must be tilted 45 degrees from the vertical axis on the projected plane 
normal to the transmission axis. In option 2, the wires of BS WG projected on the same plane as 
the projection plane for WG 2 must be perpendicular to those of WG 2. Then, the beam taken by 
WG 1 is transmitted and the beam taken by WG 2 is reflected by BS WG, and those are led to 
output WG A. The wire of output WG A must be vertical or horizontal to get interference terms 
of the two beams. 

To extract the maximum ability of MuFT, direct detectors like a bolometer, which performs 
broadband intensity measurement, must be mounted on the DeS unit. One of two bolometers 
detects transmitted signals through the output WG, while the other detects reflected signals. The 
difference in reflected and transmitted signals reproduces the full mutual coherence function since 
the phase difference in the interferograms between two signals is ir. Also, unwanted DC components 
are removed by this process. The DC fluctuations originating in the atmospheric fluctuation and 
the thermal fluctuation of the detector systems is mostly common to the transmitted and reflected 
components. Therefore, these fluctuations are supposed to be removed by taking the difference of 
each signal. This technique is also applicable to MP-FT. 

The PSI system analyzes the acquired data. Also, this unit performs mechanical control of MuFT. 
Baseline length vector source tracking and monitoring, mirror control, WG control and data acqui¬ 
sition are performed. Data analysis is also performed by the PSI system. Details of data analysis 
are explained in a forthcoming paper. 
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3.B. Detailed description of observables of MuFT 


In this section, fundamentals of MuFT are fully described by treating electro-magnetic waves as 
vector waves. Suppose the electric field of the incident electromagnetic waves emitted from sky 
position 6 on the source plane at frequency v described by E(Q,v,t) = (Ei,E 2 ) where E\ and E 2 
are orthogonal components of the electromagnetic waves. They can be written as, 


Ei 


E 2 


€i(Q, v) COs(27Tl4 + (5i) 
€2 COs (27 TVt + 82) 


\ exp[2vrjTt + jSi] + exp[-2 irjvt - j8{\ 
£i(0,v) - 2 - 

to \ exp[27 rjvt + jS 2 \ + exp[-2v xjvt - jS 2 ] 

e 2 ( 0 , v) -7T- 


( 11 ) 


The phases for the fixed frequency, hi and h 2 , can be treated as time-independent variables. In option 
1, WG1 and WG2 reflect the same polarization components of the incident electromagnetic waves. 
Therefore, option 1 measures the mutual coherence of the same component of the electromagnetic 
waves. In option 2, the components reflected by WG 1 and WG 2 are orthogonal to each other. 
Therefore, option 2 measures the mutual coherence of the orthogonal polarization components of 
the electromagnetic waves. For each option, there are two independent possible alignments of WGs, 
as explained in the previous subsection. The observed intensity transmitted through WG A for each 
option is expressed as 


/o P i(b,x) 


Iop2(t>,x) 


l/o / A ' (fl)f H 2 + eXp 


—27 xju ( b • 6 — 


2x 


+ exp 

-f 

2 Jn 


27 Tju (b-e-^pj 

f MO) 


d 1 / d 2 e, 


eiei + e 2 e 2 


+eie 2 exp[j(hi - h 2 )] exp 
+eie 2 exp [~j(5i - h 2 )] exp 


-2t rju ( b • 6 - — ^ 


27 Tjv 


dv) d 2 e, 


( 12 ) 


where time averages were performed and A U {Q) expresses a beam pattern of the observational 
system at frequency v toward 6. The factors of 1/4 and 1/2 appearing in front of each integral 
come from the transmissivity and reflectivity at BS WG and output WG. Stokes Parameters are 
defined in the usual manner 


1(0, v) = [ei(6,u)e 1 (6,u) + e 2 (0,u)e 2 (0,u)], 

Q(0, v) = [ei(0,v)e 1 (0,u) - e 2 (0,v)e 2 (0,v )\, 

U(0,v) = 2 [e\(0, v)e 2 (0, v) cos h] , 

V(0, v) = 2 [ei(0, zz)e 2 (0, v) sin h], (13) 


where 8 = 5 2 — hi is the phase difference between the orthogonal components of the incident 
electromagnetic waves. By rewriting Eq.(12) in Stokes parameters, we have 


Iorl = 4 in 


A v (0) 


l(0,v)±Q(0,v) 


1 + cos 27 t—( b • 0 — 2x) 
c 


du \ d 2 0, 
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( 14 ) 


lop 2 = ~ | j A u (G)^ 1 (0, v) + (u(0, v) cos[2vr^(b • 6 - 2x)\ 


±V(0, v) sin[27r— (b ■ 6 — 2x)\ 


du } d 2 e, 


where ± corresponds to the alignment of the WGs, that is to CW or CCW alignments described 
in the previous subsection. These results show that the MuFT is capable of measuring all intensity 
distributions of 4 Stokes parameters in a wide band, eg. 1(0, v), Q(0, u),U(0, v) and V(0, v). When 
the incident light is unpolarized, <5i and 62 are uncorrelated, and therefore, option 2 output is zero. 
For unpolarized sources, the intensity distribution can be measured by using option l(0hta(2001) 21 
, Hattori et al.(20 00) 22 ). Details of experimental results of polarimetry with MuFT are presented 
in a forthcoming paper. 


3.C. Spatial and spectral resolutions 

Spectral resolution is determined by the maximum internal time lag, tq, as A u ~ 1 /tq. Spatial 
resolution for a fixed wavelength, A, is defined in the same way as that of an conventional interfer¬ 
ometer, and is limited by the maximum baseline length to A /B max . Since the wavelength coverage 
is very wide in the case of MuFT, the maximum spatial resolution is determined by the spatial 
resolution for minimum observable wavelength. 

3.D. The field of view of the MuFT 

Suppose broadband point-source observations with two apertures using MuFT. Define the center 
of the field of view (FOV) as the origin of the source coordinate. When the source position is 9, 
there is a geometrical time lag of r g = b • 9/c. For simplicity, we assume that 6 is parallel to b and 
b • 6 = b6. As explained in Sec.2.B, the central burst of the interferogram for this source appears 
at r = — T g . To detect the central burst for this source, to must be larger than r g . This condition 
provides a limit on the FOV of MuFT of 9 < (cro)/6. A similar constraint on the FOV is deduced 
in a different way. As discussed in the previous subsection, the finiteness of the maximum internal 
time lag leads to the spectral uncertainty of Av ~ 1/ro- The frequency uncertainty results in a 
phase uncertainty of 5(f = 2nAi , b9/c. To be able to observe a visibility function without significant 
coherent loss, 5f> should be less than n/2. This condition provides almost the same limit on the 
FOV as shown in the above inequality. The limit shown above constrains the FOV only in the 
direction parallel to the baseline vector. In the direction perpendicular to the baseline vector, the 
above condition does not constrain the FOV. However, since the baseline vector must be rotated 
relative to the source plane for image synthesis, the same coverage of the FOV in both directions 
might be convenient. Therefore, the maximum available pixel number of the detector array installed 
on the focal plane of the MuFT is provided by (cTo/bO p ) 2 , where 9 P is the FOV of one detector 
pixel. 

The analysis of the data measured by an off-axis pixel is not self-evident. Since the position of 
the central burst is shifted from the scanning center, the Fourier transformation of the obtained 
interferogram in t from —tq to tq causes a phase shift of 27 TVT g /c. Therefore, the obtained phase 
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should be shifted back by ‘I'nvTg/c. There are other non-trivial problems related to using the multi¬ 
pixel detector array on MuFT which must be studied to optimize methods. Rinehart et al. (2004) 10 
have been studying the application of the multi-pixel detector array on the double Fourier in 
laboratory experiments in optical wave bands using an optical CCD camera. They have successfully 
detected interferograms in each pixel which are each observing different FOV. 

3. E. Phase errors, efficiencies and noise 

There are four sources of phase errors: (1) sampling error due to shift in zero optical path difference; 
(2) uncertainty of (u, v ) points due to finiteness of spectral resolution; (3) uncertainty of (u, v ) points 
due to pointing error. (4) Atmospheric phase fluctuation. The first one can be compensated by zero- 
path calibration of the interferometer. The second one provides a limit on the FOV, as discussed 
in the previous subsection. The third one can be corrected by pointing observations of compact 
sources with a certain interval. Atmospheric phase error can be corrected when atmospheric opacity 
and phase error are correlated. 

As seen in Eq.(12) and Eq.(14), optical efficiency of MuFT for a single output port is 1/8 for 
option 1 when unpolarized radiation is observed, combination of two output ports, such as that 
used for detecting reflection and transmission of WG A, would give an efficiency of 1/4. Detecting a 
signal traveling toward WG B would give an efficiency of 1/2. Furthermore, setting another MuFT 
under the first one to detect the radiation transmitted through WGs 1 and 2 would give an efficiency 
of 1. 

Sources of noise in observed interferograms are atmospheric radiation, radiation from optical 
components, internal detector noise and other excess noise such as readout electronics noise. Since 
most of these noises are incoherent, they do not create coherent peaks in an interferogram but 
appear as random noise. Common amplitude fluctuation of atmospheric emission could be removed 
by taking the difference of the interferograms obtained by signals transmitted through and reflected 
by output WG. Radiation from optical components causes a small amount of coherent noise signal 
in an interferogram. Most of the noise emission from the optical components which propagate back 
to the input wire grids is ejected to the sky. However, the direction of noise emission propagated 
from optical components has a finite distribution, and finite errors in alignment of wire grids 
and optical components are not escapable. Therefore, a small amount of radiation returns to the 
detector. This results in an autocorrelation noise emission signal from the optical components in 
the interferogram. This can be removed by subtracting the background interferogram obtained by 
observations of the sky without source. However, it is impossible to observe exactly the same sky 
as during source observations without the source. We have to adopt data obtained by observations 
of the sky without the source as the background data. 

4. Advantages of MuFT 

One of the most prominent advantage of MuFT is that one can use direct detectors like a bolometer 
as an interferometer focal plane detector. This give MuFT the following three significant advantages. 
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High sensitivity and wide bandwidth interferometer in FIR As described in the introduc¬ 
tion, the sensitivity of direct detectors is higher than that of heterodyne receivers in the frequency 
range of v > 100 GHz, that is A < 3 mm. Therefore, the sensitivity of MuFT could be higher than 
an interferometer based on heterodyne receivers, especially in submillimeter and FIR wave bands. 
Performing wide band observations with direct detectors is easily achieved. For example, bolometers 
which have a sensitivity in extremely broadband have been used as focal plane detectors of FTS in 
FIR wave bands. Since the WGs have uniform reflectivity and transmissivity for all frequencies in 
millimeter, submillimeter and FIR, the MuFT is suited to broadband observations. 

Wide field of view Constructing large format heterodyne focal plane arrays is not yet feasible, 
since local oscillators and associated electrical circuits for each detector are bulky and not amenable 
to automated production using integrated circuit technology. Therefore, not a single detector is used 
in radio interferometers, including ALMA. This severely limits the FOV of the interferometer. On 
the other hand, large format focal plane arrays of direct detectors are already operating. New types 
of direct detectors aiming at a dramatic increase in the number of pixels are lying developed. 5,6 
Since MuFT can use direct detectors as interferometer focal plane detectors, it could extend the 
FOV of the millimeter and submillimeter interferometers dramatically. Details are described in 
Sec.3.D. 

Large dynamic range When source intensity distribution is independent of frequency, we can 
sum up all images obtained at different frequencies into a single image, as described in Sec.2.A. In 
other words, a single baseline observation is equivalent to observations with many combinations of 
baseline interval in the case of a single frequency interferometer. Long wavelengths observe overall 
source structures with low spatial resolution. Short wavelengths observe fine structures with high 
spatial resolution. A large dynamic range of spatial resolution is covered even by a single baseline 
observation. A similar technique has been addressed as MFS for heterodyne interferometer(Conway 
et al.(1990) 14 ). This is a technique for combining maps obtained for different frequencies by different 
interferometers into a single image. 

Linear and Circular Polarimetry The MuFT can measure all 4 components of Stokes pa¬ 
rameters (Hattori et al.(2000) 22 ), as explained in sec. 3.B. This is because MuFT uses WGs as 
beam-splitter. Since a double Fourier uses dielectric film as a beam-splitter, it cannot perform 
polarimetry. 

5. Summary 

We proposed a new type of bolometric interferometer named MuFT by applying MP-FT to aper¬ 
ture synthesis. Fundamentals of imaging, spectroscopy and polarimetry with this instrument were 
developed. The MuFT is a system which permits imaging and spectroscopy in a wide band by 
combining the Wiener-Khinchine Formula 13 which makes it possible to extract a spectrum from 
the auto-correlation and van Cittert Zernike Formula 13 which allows imaging from the mutual 
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correlation. The concrete composition of this equipment was proposed. By combining wire grid 
beam-splitters adequately, source intensity distributions of four Stokes parameters can be acquired 
in a wide band. Fundamental restrictions in practical use were discussed. 

We are planning to perform the observations from good mm and submillimeter observation sites, 
such as the Nobeyama Radio Observatory in Japan, the Atacama desert in Chile, South Pole etc.. 
The possibility of having a super-wide FOV by mounting the focal plane array of bolometers, SIS 
photon detectors and transition edge sensors is attractive for future applications. Optical designs 
of MuFT for using detector arrays is going to be advanced. 

Applying this technique to a space-borne mission is one of the best possibilities for extracting 
the maximum ability of MuFT, since there is no restriction on the bandwidth from atmospheric 
absorptions. Mather et al. 23 has been proposing a space-borne FIR observatory based on this kind 
of technique. The future application of this technique to observations from space could open new 
and interesting possibilities in FIR astronomy. 
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List of Figure Captions 

Fig. 1.Coordinate systems and designations used in observation of extended source. Source is as¬ 
sumed to be very far away. 

Fig. 2. Simplified schematic diagram of MuFT. Wave-front division of source light is performed 
by LiC. These are combined by FI, through optical systems. Interferogram is measured by DeS. 

Fig. 3. Example of FI unit. WG position is Optionl. Thin dashed lines show light paths. Arrows 
attached to the dashed lines indicate the propagation direction of the light beams. 

Fig. 4. Two basic configurations of wire grids 1 and 2. The left figure is Option.1 (WGs are 
parallel) and the right figure is Option.2 (WGs are perpendicular). Gray lines are light paths. 
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Fig. 1. Coordinate systems and designations used in observation of extended source. 
Source is assumed to be very far away. 
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Fig. 2. Simplified schematic diagram of MuFT. Wave-front division of source light is 
performed by LiC. These are combined by FI, through optical systems. Interferogram 
is measured by DeS. 
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Fig. 3. Example of FI unit. WG position is Optionl. Thin dashed lines show light 
paths. Arrows attached to the dashed lines indicate the propagation direction of the 
light beams. 



Fig. 4. Two basic configurations of wire grids 1 and 2. The left figure is Option. 1 
(WGs are parallel) and the right figure is Option.2 (WGs are perpendicular). Gray 
lines are light paths. 
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